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Motivation 

Rapid  Communication  is  essential  in  today's 
world. 

Understand  the  dynamics  of  flow  by  creating  a 
communication  network  model. 

Modeling  flow  on  a  communication  network  will 
allow  us  to: 

—  Describe  normal  and  congested  flow  on  large 
communication  networks. 

-  Predict  changes  in  flow  pattern  due  to  changes  in 
the  spatial  density  and  per-link  traffic. 
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Communication  Networks 

•  A  communication  network  is  a  global  system  of 
interconnected  networks,  both  big  and  small. 

•  Packet  switching  network  because  all  data  traffic 
is  broken  down  into  data  chunks  called  packets. 

•  Everything  traveling  on  a  communication  network 
is  called  a  packet. 
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Discrete  Communication  Network 


•  Use  graph  theory  to  describe  the  connectivity  of  a 
network,  where  a  graph  is  composed  of  nodes 
and  links. 

•  Information  travels  along  links  connecting  nodes. 

•  The  graph  is  undirected;  information  travels  in 
both  directions. 

•  The  nodes  are  routers  in  this  model 

•  Routers  will  act  as  both  host  and  switch 
computers. 

•  Host  computers  is  where  information  enters 
(source),  and  exits  (destination)  a  network. 
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Route  Matrix 


The  Route  matrix  depicts  the  global  state  of  a 
network. 

Describes  how  to  direct  packets  from  source  to 
destination. 

Each  entry  describes  the  next  appropriate  router 
a  packet  will  take  along  its  path. 

Routes  are  pre-determined  by  an  optimal  path 
algorithm. 
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Queue  Dynamics 

•  Each  router  contains  a  queue  of  buffered 
packets. 

•  FIFO  unlimited  memory  buffer 

—  Packets  enter  queue  in  one  of  two  ways 

-  Packets  flowed  from  another  router  along  a 
connecting  link. 

-  Generated  at  a  router  according  to  some  packet 
generator  rate  appropriate  to  the  router. 

•  Upon  being  generated,  a  packet  is  given  a 
destination  by  its  originator. 

•  Once  a  packet  reaches  its  destination  it  is 
delivered  and  exits  immediately. 
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Flow  equation  (Arrivals  Age  zero) 


oo  N 

b(j ,  d,  0,  t  + 1)  =  ^  ^  8jRJ3{i ,  d,  a ,  r)  +  v(j,  d,  t) 

a= 0  i— 1 
j*d 

b(j,d,a,T),  number  of  buffered  packets  at  node  j, 
with  destination  J,  age  a  at  time  t 

(5{i,  d,  a,  z)  is  the  sending  rate  of  packets  being  sent  from  router 
i  to  router  j  toward  its  destination  d  at  time  r 

v(j,  d,  t )  is  the  number  of  new  packets  entering  the  network  at 
node  j 
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Evolution  Equation 

•  Describes  how  packets  are  aging  in  the  buffer. 

b(j ,  d,  a  + 1,  t  + 1)  =  b(j ,  d ,  a ,  r)  -  fl(j\  d ,  <2,  r) 

/?(y ,  a,  r)  determines  the  number  of  packets  that  are  sent  outward 
from  node  j  with  destination  d 
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Discrete  Conservation  of  Packet 

Equation 

The  number  of  buffered  messages  at  router  j  with 
destination  d  at  time  t 

oo 

n(J,d,T)  =  '^b(J,d,a,T) 

<2=0 


N  N 

n(j,  d,T  +  l)-n(j,d,T)  =  ^  b(j,  d,  a,  r  + 1)  -  ^  b(j,  d,a,r ) 

a= 0  a=0 


Network  Science  Center 

West  Point 


Discrete  Conservation  of  Packet 

Equations 


OO  00  N 

n(j,  d,T  + 1)  -  n(j,  d,  t )  =  J3(j,  d,a,T)  +  ^^  diRij3(i,  d,  a,  t) 

<2=0  <2=0  i= 1 


+  v(j\  d,r)~Yj  P(j,  j,  a,  r) 

<2=0 
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Continuum  Network  Model 

•  To  view  this  model  as  a  flow  model,  well  discuss 
the  collection  of  routers  as  opposed  to  one.  The 
collection  of  routers  create  a  Voronoi  Diagram. 
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Voronoi  Diagram 

The  spatial  location  of  each  router  is  of 
importance. 

Each  router  serves  a  particular  coverage  area  of 
users  who  are  sending  and  receiving  packets  in  an 
area 

Associate  each  router  with  a  physical  spatial 
location 

Each  router  shares  physical  boundaries  with 
another  set  of  routers  so  that  packets  moving 
through  the  network  will  pass  through  physical 
boundaries  around  routers. 

Network  Science  Center 

West  Point  ^ 


Voronoi  Diagram 


x ,  =  Location  of  router  / 

J  J 

y  =  Location  of  destination  d 


boundary  of 
x(j) - ^x(k) 


Each  Xj  generates  a  Voronoi  polygon  V (x  . ),  a  tessellation  in  the 
plane.  The  Voronoi  polygon  V (jc  . )  is  the  set  of  points  closer  to 
the  point  (router)  x. ,  than  to  all  other  routers  in  the  diagram.  A 
collection  of  Voronoi  polygons  is  called  a  Voronoi  Diagram. 
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Continuum  Description 


Let  Vp  be  a  Voronoi  Diagram  and  dVp  be  its  boundary 
The  density  of  packets  for  Voronoi  polygon  V (x  ) 


p(Xj,y,t) 


njj^d^T) 

v(xj) 


The  number  of  packets  buffered  in  the  Vorono  Diagram  Vp 

j^p 
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Continuum  Description  Cont. 


The  evolution  equation  for  the  density  of  packets  in  a  Voronoi  Diagram 


L 


p(x,  y, t  +  At)  -  pjx,  y,  t)  ^  =  y  n(j,  d,  r  + 1)  -  n(j,  d, r) 


At 


JeP 


At 


Take  the  limit  as  At  goes  to  zero,  first  term  becomes 
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Flux 


The  flow  vector  from  router  x,  to  router  x,  is 


j 


x .  -  X  . 

1  J 


x .  -  X  . 

L  J 


fi(j,d,a,r ) 


The  outflow  through  boundary  element  <3F;  m 


- 1  x,.  -  xy 


"77” - -P(j’d,a,T)-nlm  =®°(xl  m,y,t)-nl  m  dVlm 


At 


x .  -  x . 

^  J 


Similar  derivation  fo  the  incoming  flow  terms 


oo 


a=0  jeP 
i<£P 


1  Xj  —  X  ■  i 

Su,mt,m)5j,R,.d— - p(i,d,a,T)-nlm  =0  (xlm,y,t)-n 


At 


x .  -  X  • 

1  J 


l,m 


dV, 


l,m 
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Flux  Cont. 


The  total  flux  of  packets  entering  and  exiting  a  boundary  dVt  m 


> TO ' ni,m \dVhm I  =  3)° (*,,m ,y,t)- nlm  dVlm  -<t>‘ (xhm,y,t)-n,JdVhm 


Total  flux  of  packets  in  Voronoi  Diagram  Vp  with  destination  y 

-n,m  8V,m 


/  j  V  /, m 

l<=dP 
m£dP 


Written  in  the  continuum  limit  and  using  the  divergence  theorem 
the  total  flux  becomes 


Z0(Xm’->v) 


n 


Lm 


dV, 


l ,  m 


l&dP 

m<£dP 


\gv  <&{x,y,t)-nds 
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Source  and  Sink 


The  source  and  sink  in  their  continuum  limit 


oo 


a=0  jeP 


1  r 

p(j,d,a,r)  — >•  cr(x,t)dV 

JVn 


At 
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Continuity  Equation 


Putting  the  continuum  limits  together  we  have  the  conservation  of 
packets  equation  in  a  Voronoi  Diagram  Vp 


dp 

—  +  V-®(x,y,t)  =  y(x,  y,  t )  -  cr(x ,  t) 
ot 
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One  Dimensional  Network  Flow 

Model 

One  dimensional  Network  flow  model  with  x=0 
and  a  destination  x=y. 

Analyze  inner  nodes 

Continuity  equation  in  one  dimensions 

dt  dx 


p(x,  y,  t)  if  p(x ,  y,  t )  <  <Dmax  (x,  y,  t ) 
(x,  y,  t )  if  p(x,  3/,  0  >  <Dmax  (x,  y,  t ) 


max 
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Non-Saturated  Flow 


t - 1 - 1 - 1 - 1 - 1 - 1 - r 


Flow  Movement  In  Time 
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End  of  Flow 
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Example:  Flow  with  Saturation 


Initial  Condition 
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End  of  Flow 
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Interruptions 

•  Disturbance  in  the  flow 

-  Limited  bandwidth,  link  capacity  drops  to  a  lower  value 

-  Router  (grid  point)  is  down  for  some  time. 

•  Example 

-  Source  node  has  constant  flux  of  packets  below  link 
capacity. 

-  One  of  the  grid  points  has  limited  bandwidth  for  a  while. 

-  One  of  the  inner  grid  points  has  a  source  term. 
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Interruptions 


t=0 

t=10 

t=30 


30  40  50  60  70  80 
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End  of  Flow 


t=80 
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Current  Work 

Currently  route  matrix  is  static. 

Make  route  matrix  dynamic  at  every  node 
(changes  because  of  link  weight,  upstream  traffic, 
etc. 


J?  —  /  ^ -  Before 

y'W  “ 1 

Rj,d  (o  *  Rj,d  (t + 1)*  n°w 
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Current  Work 


Also  leads  to  probabilities  in  which  nodes  packets 
will  take  next. 


Probability  of  going  to  node  /  to  node  j  with 


destination  d. 

P(j,d,t ) 


Destination 
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Questions!!!!!! !!!!!! 


Network  Science  Center 

West  Point  '2$ 


